Translational Entanglement of Dipole-Dipole Interacting Atoms in Optical Lattices 
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We propose and investigate a realization of the position- and momentum-correlated Einstein- 
Podolsky-Rosen (EPR) states [Phys. Rev. 47, 777 (1935)] that have hitherto eluded detection. The 
realization involves atom pairs that are confined to adjacent sites of two mutually shifted optical 
lattices and are entangled via laser-induced dipole-dipole interactions. The EPR "paradox" with 
translational variables is then modified by lattice-diffraction effects, and can be verified to a high 
degree of accuracy in this scheme. 

PACS numbers: PACS numbers: 03.65.Ud, 34.50.Rk, 34.10.+X, 33.80.-b 



The ideal Einstein-Podolsky-Rosen (EPR) Q] state of 
two particles — 1 and 2, is, respectively, represented in 
their coordinates or momenta (in one dimension), as fol- 
lows, 
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The "paradox" is in the fact that given the measured 
values of x\ or p\ of particle 1, one can predict the 
measurement result of x 2 or p 2 , respectively, with ar- 
bitrary precision, unlimited by the Heisenberg relation 
Ax 2 Ap 2 > Ti/2. In other words, the ideal EPR state is 
fully entangled in the continuous translational variables 
of the two particles. Approximate versions of this transla- 
tional EPR state, wherein the J-function correlations are 
replaced by finite- width distributions, have been shown 
to characterize the quadratures of the two optical-field 
outputs of parametric downconversion 0, and allow 
for optical continuous- variable teleportation Q. More 
recently, translational EPR correlations have been an- 
alyzed between dissociation fragments of homonuclear 
diatoms whereas interacting atoms in Bose-Einstein 
condensates have been shown to possess translational- 
internal correlations .6]. Yet the fact remains that the 
original EPR state has eluded detection for nearly 70 
years. Our goal is twofold: (i) propose an experimentally 
feasible scheme for the creation of translational EPR cor- 
relations between cold atoms that are confined in opti- 
cal lattices @ and coupled by laser- induced dipole-dipole 
interactions (LIDDI) 0, [jj, QiJ ; (ii) study the qualitative 
modifications of such correlations due to particle diffrac- 
tion in lattices, which have been hitherto unexplained. 
The LIDDI has been proposed as a means of two-atom 
entanglement via their internal states, for quantum logic 
applications The ability of LIDDI to influence the 
spatial and momentum distributions of cold atoms in cav- 
ities traps and condensates 01, has been investi- 
gated extensively. 

To realize and measure the EPR translational correla- 
tions of material particles, one must be able to accom- 
plish several challenging tasks: (a) switch on and off the 
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FIG. 1: Proposed scheme of the translational EPR: two over- 
lapping optical lattices displaced from each other in the y 
direction by I, are sparsely occupied by two kinds of atoms. 
Each of the two kinds of atoms feels a different lattice; the 
shaded regions depict the energy minima (potential wells) of 
the lattices. 



entangling interaction; (b) confine their motion to sin- 
gle dimension, and (c) infer and verify the dynamical 
variables of particle 2 at the time of measurement of par- 
ticle 1. The latter requirement is particularly hard for 
free particles, since by the time we complete the predic- 
tion for particle 2, its position will have changed. In [B| 
we suggested to overcome these hurdles by transforming 
the wavefunction of a flying (ionized) atom by an elec- 
trostatic/magnetic lens onto the image plane, where its 
position corresponds to what it was at the time of the di- 
atom dissociation. In this Letter we propose a different 
solution: (a) controlling the diatom formation and disso- 
ciation by switching on and off the LIDDI; (b) controlling 
the motion and effective masses of the atoms and the di- 
atom by changing the intensities of the lattice fields. 

System specification: Let us assume two overlapping 
optical lattices with the same lattice constant a, as in 
Fig. □ The lattices are very sparsely occupied by two 
kinds of atoms, each kind interacting with only one of 
the two lattices. This can be realized, e.g., by assum- 
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ing two different internal (say, hyperfine) states of the 
atoms llj. For both lattices, the y and z directions are 
very strongly confining (realized by strong laser fields), 
whereas in the x direction the lattice can be varied from 
moderately to weakly confining. Thus, the motion of 
each particle is confined to the x direction. For each 
direction we assume that only the lowest vibrational en- 
ergy band is occupied. Initially, the potential minima of 
the lattices are displaced from each other by an amount 
I <C a in the y direction. This enables us to couple the 
atoms of the two lattices, along y using an auxiliary laser 
to induce the LIDDI. We assume the auxiliary laser to 
be a linearly polarized traveling wave with wavelength 
Ac, moderately detuned from an atomic transition that 
differs from the one used to trap the atoms in the lattice. 
The auxiliary laser propagates in the x direction and its 
electric field is polarized in the y direction. The LIDDI 
potential for two identical atoms has the form y| 
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Here Vc = a 2 k 3 I c / {^ire^c) , where k = 2ir/\ c , I c is the 
coupling laser intensity, and the atomic dynamic polar- 
izability is a = 2uja\p\ 2 /[Ti(lj a — u> )], /i being the dipole 
moment element, loa the atomic transition frequency, 
and u> — kc. The position-dependent part F§(kR) is a 
function of R, the distance between the atoms, and 9, 
the angle between the interatomic axis and the wavevec- 
tor of the coupling laser. Since I <C 2a, V dd (R) has a 
pronounced minimum for atoms located at the nearest 
sites, R~l, where V dd {R) ~ -Vc(Ac/0 3 /(4tt 3 ). Under 
these assumptions, we can treat the system as consist- 
ing of pairs of "tubes" , that are oriented along x, either 
empty or occupied. Only atoms within adjacent tubes 
are appreciably attracted to each other along y, due to 
the LIDDI. 

EPR states: We now focus on the subensemble of tube- 
pairs in which each tube is occupied by exactly one atom. 
In the absence of LIDDI, the state of each atom can be 
described in terms of the Wannier functions \xj) EH that 
are localized at lattice sites with index j and may hop to 
the neighboring site at the rate \\\ op /h, where Vhop = 
(Xj\H\ a t\xj+i) i an d -Hiat is the lattice Hamiltonian. In 
a ID lattice, H\ &t = (U /2) cos (2irx/a) + p 2 x /(2m), m 
being the atomic mass. The hopping rate is related to 
the energy bandwidth of the lowest lattice band Vb by 
Vb ~ 4 1 Vhop | (for exact expressions see 0]). For a 
shallow lattice potential (Uo < 15E ICC ) we may use the 
approximate formula Vhop ~ E roc exp(— 0.26 Uo/E iec ), 
where the recoil energy is 

Let us switch on the LIDDI, so that |Vh op | -C |Vdd|- 
Then the ground state of such a tightly bound diatom can 

be approximated by \ipo) oc J2j \Xj^)\Xj^)- This means 
that when particle 1 is found at the jth site of lattice 
1, then particle 2 is found at the jth site of lattice 2, 
with position dispersion given by the half-width a of the 




FIG. 2: Joint probability distribution of the positions of 
two lithium atoms in adjacent optical lattices, prepared in a 
diatom state as specified in the text, using the ground state of 
the external harmonic potential with half- width of ge — 6a 
and temperature of 10 nK. Inset: Position probability of atom 
2 in the state above, conditional on atom I being measured 
at site (full line). Dashed line: Gaussian approximation of 
the Wannier function with the half- width a = 0.14a. 



atomic (Gaussian-like) Wannier function in the lowest 
band, a 2 = hX^ / (Airy/mUo) . To next order in Vh op /Vdd, 
the nonzero probability of atoms to be located at more 
distant sites changes the diatomic position (separation) 
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The states of the tightly bound diatom form a separate 
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band whose bandwidth is V L 
lowest atomic vibrational band 
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The diatomic hopping 
potential V^p 1 ' can be found by assuming that the two 
atoms consecutively hop to their neighboring sites, i.e., 
the state change U^Hu^x _ k.M \u.( 2 ) 

either via \x^ 
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via \ X f)\xf ] ) - \xf)\xfl) - IxSiMxf+T)- By adia- 
batic elimination of the higher-energy intermediate states 

one obtains V h ( D 2 p at) « 2V 2 op /V dd . 

To realize a momentum anti-correlated EPR state, the 
temperature of the system must satisfy fc^T <C Vg . 
The dependence of the momentum anti-correlation on 
temperature is given by Ap^_/(2m[.g at ') s» |fcgT, where 
we have introduced the sum-momentum spread Ap + and 
(analogously to atomic effective mass |14|) the two-atom 
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We then ob- 



tain Ap 2 ^ 



n 2 \v dt 



rfcsT. 



Although the values Ax 2 _ and Ap 2 ^ estimated above 
refer to the respective peak widths, they principally dif- 
fer from the position and momentum uncertainties of 
free particles: due to the lattice periodicity, the position 
and momentum distributions have generally a multi-peak 
structure. The two-particle joint position distribution of 
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FIG. 3: Conditional probability of the momentum of atom 2 
after the momentum of atom 1 has been measured for lithium 
diatoms prepared as in the text (the measured value pi = piM 
is indicated with an arrow). The dashed line corresponds 
to the marginal probability distribution of momentum p2 ir- 
respective of the momentum of atom 1 at the temperature 
T — 100 nK. The half-width of each peak is equal to 1/s of 
Eq. (|1J. Inset: joint probability distribution of the atomic 
momenta in the state above with T = 100 nK. 



the ground state is a chain of peaks of half-width a sep- 
arated by a; the peaks are located along the line x 2 = x\ 
(Fig. [21 . The corresponding joint momentum distribu- 
tion spreads over an area of half- width h/(2a) and con- 
sists of ridges in the direction p 2 = —pi- These ridges 
are separated by 2-Kn/a, and for a lattice of N sites, the 
half-width of each ridge is ttT%/{No) (Fig. [^J. 

To evaluate how "strong" the EPR effect is, we com- 
pare the product of the half-widths of the position and 
momentum peaks in the tightly bound diatom state de- 
scribed above with the limit of the Heisenberg uncer- 
tainty relations, AxAp > h/2, defining the parameter s 
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A value of s higher than 1 indicates the occurrence of 
the EPR effect; the higher the value of s, the stronger 
the effect. Strictly speaking, because of the multi-peak 
momentum distribution, one should not use the original 
form of Heisenberg uncertainty relations but a more gen- 
eral relation, as discussed, e.g., in 0, that distinguishes 
the uncertainty of a few narrow peaks from that of a 
single broad peak. However, even the simple half-width 
of the peaks is a useful measure of the EPR effect. In 
order to maximize s, we must adhere to the trade-off be- 
tween making Ax- as small as possible, by decreasing 
l^hop/Kidl) and making Ap + as small as possible, by in- 
creasing |Vho P /Vdd|- The optimum value of s generally 
depends on the lowest available temperature of the di- 
atom, as detailed below. 

EPR state preparation: Cooling down the diatomic 
system to prepare the EPR state is a non-trivial task. 



We suggest to attack the problem by a three-step ap- 
proach, (i) Let us switch off both the LIDDI and the 
x-lattices, switch on an external, shallow, harmonic po- 
tential in the x direction, and cool the x-motion of the 
atoms down to the ground state of the external poten- 
tial. The width erg of the ground state should be several 
times the lattice constant; it is related to the desired mo- 
mentum anti-correlation by o~e ~ ^/(V2Ap + ). The tem- 
perature must be T <§C h 2 / '{Araks &%)■ (ii) A weak lattice 
potential in the ^-direction is then slowly switched on, so 



that the state becomes 
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where the co- 



efficients Oj ~ exp[— (j — j ) 2 a 2 /(4a E )] are Gaussians 
localized around the minimum of the external potential, 
(iii) We switch on the LIDDI and change the sign of the 
external potential, from attractive to repulsive, acting to 
remove the particles from the lattice. The two parts of 
the wavefunction would behave in different ways. The 
paired atoms, corresponding to the part of the wavefunc- 
tion V . a 2 |Xj^ ) IXj 2 ^ ) > move slowly because of their large 
effective mass m^ at ' ) , whereas single (unpaired) atoms, 



"off 

because of their smaller effective mass, 
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m cff < m; ff - 
|Vdd|/(2Vh p)TO e ff, are ejected out of the lattice and sep- 
arated from the diatoms as glumes from grains. The 
paired atoms remaining in the lattice are then in the state 

~ exp[-(j - jo) 2 a 2 / {2a 2 ,)]\xf ) )\xf ') wherein positions 
are correlated with uncertainty Ax+ « and mo- 

mentum uncertainty Ap + w h/ Ax + . At higher tempera- 
tures the atoms are not cooled to the ground state of the 
external potential and the momentum anti-correlation 
becomes Ap + « %/{ ^2a E t&nh[h 2 /(2o- 2 E mk B T)}}. The 
parameter s of Eq. (0J can then be estimated as 



V2a 



tanh 
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This equation enables us to select the optimum external 
harmonic potential (specified here by o~e) such that the 
parameter s is maximized, under the constraint of the 
lowest achievable temperature T. 

The small effective mass of unpaired atoms allows us to 
cool them individually, restricting their cooling to tem- 
peratures higher than that corresponding to the bottom 
of the diatomic band. The price is, however, that most 
of the atoms arc discarded and only a small fraction re- 
mains in the diatom state. Specifically, out of the total 
number of tube pairs occupied by two atoms, a fraction 
of ~ a/o~ E will remain in the bound diatom state. The 
different behavior of the paired vs. unpaired atoms in a 
periodic potential is a sparse-lattice analogy of the Mott- 
insulator vs. superfluid state of the fully occupied lattice, 
recently observed in Ref. [l7j . 

Measurements: After preparing the system in the EPR 
state, one can test its properties experimentally. To this 
end we may increase the lattice potential Uq, switch off 
the field inducing the LIDDI, and separate the two lat- 
tices by changing the laser-beam angles. By increasing 



4 



Uq, the atoms lose their hopping ability and their quan- 
tum state is "frozen" with a large effective mass: the 
bandwidth Vg decreases exponentially with Uq and the 
effective mass increases exponentially, so that the atoms 
become too "heavy" to move. One has then enough time 
to perform measurements on each of them. 

The atomic position can be measured by detecting its 
resonance fluorescence. After finding the site occupied 
by atom 1, one can infer the position of atom 2. If this 
inference is confirmed in a large ensemble of measure- 
ments, it would suggest that there is an "element of real- 
ity" [l| corresponding to the position of particle 2. The 
atomic momentum can be measured by switching off the 
x-lattice potential of the measured atom (thus bringing 
it back to its "normal" mass m): the distance traversed 
by the atom during a fixed time is proportional to its mo- 
mentum. One can test the EPR correlations between the 
atomic distributions occupying the two lattices: a large 
number of pairs would be tested in a single run. The 
correlations in x and anti-correlations in p would be ob- 
served by matching the distribution histograms measured 
on particles from the two lattices. 

Example: We consider two lithium atoms in two lat- 
tices with Al = 323 nm (corresponding to the transition 
2s-3p) and a dipole-dipole coupling field of Ac = 670.8 
nm (transition 2s-2p). The field intensities are /l = 
0.35 W/cm 2 and I c = 0.1 W/cm 2 , and the field detun- 
ings are <5l = 507l, <5c = 1007c, the decay rates being 
7L = 1.2 x 10 6 s" 1 , and 7c = 3.7 x 10 7 s" 1 . The two 
lattices are displaced by Z = 40 nm. From these values 
we get the lattice potential Uq — 7A2E rcc , the dipole- 
dipole potential of the nearest atoms Vdd = — 2.16£ , roc , 
and the hopping potential Vj lop = — 0.0355E ICC . The two- 
particle hopping potential is then ~ — 0.0012i? rec 
and the ratio of effective masses of a diatom and a of sin- 
gle atom is rn^^ /m e g ~ 30. The position uncertainty 
of atom 2, after position measurement of atom 1, is then 



Ax_ « a = 0.136a = 22 nm (see Fig. EJ). The correlated 
pairs are prepared by first cooling independent atoms 
in an external harmonic potential with the ground-state 
half- width of erg = 6a (frequency of 1.2 kHz ~ 33 nK). 
After the unpaired atoms are removed from the lattice, 
we calculate the momentum distribution for two different 
temperatures, 10 nK and 100 nK. The conditional prob- 
ability of momentum pi of particle 2, provided that the 
momentum of particle 1 was measured as p\u is plotted 
in Fig. |3 The resulting half- widths of the peaks can be 
used to find the parameter s; we have s ~ 30 for T — 
10 nK, and s m 11 for T = 100 nK. Note that in current 
optical experiments Q s < 4. 

To sum up, the proposed scheme is based on the adap- 
tation of existing techniques (optical trapping, cooling, 
controlled dipole-dipole interaction) to the needs of atom- 
atom translational entanglement. The most important 
feature of the scheme is the manipulation of the effec- 
tive mass, both for the EPR-pair preparation (by sepa- 
rating the "light" unpaired atoms from the "heavy" di- 
atoms) and for their detection (by "freezing" the atoms 
in their initial state so that their EPR correlations are 
preserved long enough). This scheme has the capacity 
of demonstrating the original EPR effect for positions 
and momenta, as discussed in the classic paper Q. A 
novel element of the present scheme is the extension of 
the EPR correlations to account for lattice-diffraction ef- 
fects. Applications of this approach to matter teleporta- 
tion [al and quantum computation with continuous vari- 
ables 18] can be envisioned. The fact that our system 
represents a blend of continuous and discrete variables 
may be utilized for quantum information-processing (to 
be discussed elsewhere). 
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